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We present a method based on three-particle azimuthal correlation cumulants for the 
study of the interaction of jets with the medium produced in heavy ion collisions at RHIC 
and LHC where jets cannot be reconstructed on an event-by-event basis with conven- 
tional jet finding algorithms. The method is specifically designed to distinguish a range 
of jet interaction mechanisms such as Mach cone emission, gluon Cerenkov emission, 
jet scattering, and jet broadening. We discuss how anisotropic flow background conti- 
butions of second order (e.g. r>|) are suppressed in three-particle azimuthal correlation 
cumulants, and discussed specific model representations of di-jets, away-side scattering, 
and Mach cone emission. 



1. Introduction 

The observation of a dip at 180° in flow subtracted two-particle azimuthal correlations 
observed in Au + Au collisions d has recently generated quite a bit of interest. Stoecker 
2 and independently Shuryak et al. El suggested the observed dip might be an indica- 
tor of the production of away-side parton induced wake field or Mach cone, with the 
width of the cone determined by the sound velocity in the produced medium. The dip 
might however also result from other processes such as large angle gluon radiation^ 
jets deflection by radial flow, or Cerenkov gluon radiation 1^1. While discrimination of 
these production mechanisms is not possible with two-particle correlations, it might be 
achieved with three particle correlations. Preliminary such three particle analyses pre- 
sented at QM05 and elsewhere by STAR El, and PHENIX0 arise as rather challenging 
tasks. Observation of three-particle correlations require large datasets. One must also 
account for the fact that measured three-particle densities consist of a superposition of 
correlated three particle signals, and combinatorial terms involving only two correlated 
or three uncorrelated particles. Finally, the interpretation of the data is also complicated 
by finite anisotropic flow, and momentum conservation effects El. 

We present and discuss the merits of an analysis technique based on cumulants. 
The definition of cumulants, measurement method, and key properties are discussed 
in Section 2. Section 3 includes a brief discussion of correlation shapes expected from 
simple models of in-vacuum jets, deflected jets, conical emission, and jet-flow. 

Some of these topics were already discussed i n El A lso note that alternative analysis 
techniques were discussed in the recent literature I ^ I IO I 
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2. Cumulant Definition and Properties 

The cumulant method makes no assumption about the background or shape of the 
signal. Cumulants are defined as statistical measures of the degree of correlation between 
measured particles. For three particles, the definition reads: 

C3(<Pi,<Pj,<Pk) = Pa{<PU<Pi,<Pk) - P2{'Pi,<P})Pl(.<Pk) ~ P2(<Pi,<Pk)Pl(tPi) 
-P2 {<Pj > fk)pi (<«) + 2/01 (<Pi)pi (<Pj)pl {fk) 

where p3(^Pi, <fj , fk) = dN/dipidtpjdip k , p 2 (tpi,ipj) = dN/dtfidfpj, and pi((pi) = dN/dipi, 
are respectively three, two, and single particle densities (normalized per event), measured 
for selected particles identified with labels i, j, and k. The definition holds whether one 
considers identical particles, or different species, or whether the integrated kinematic 
ranges are the same, overlap, or are different for the three particles measured. 

By construction, cumulants are non-positive definite quantities, and indicate the 
degree to which particles are correlated relative to Poisson processes: positive values 
indicate regions of phase space where particles are likely to be founded together, whereas 
negative values signal regions where particles are unlikely compared to uncorrelated 
particle production (Poisson). It is straightforward to show that the cumulant of a sum of 
independent processes is equal to the sum of the cumulant of each the processes thereby 
enabling modeling of particle production in terms of separate components, such as jets, 
and background particles. Note however that such a separation in a finite size system 
leads to sampling biases: cumulants are non-zero even in the absence of correlations. 

The definition Eq. 1 is of limited interest for the study of the physical processes 
involved in the particle production. More relevant is a formulation in terms of relative 
angles. With three particles, there are two such independent angles we here choose as 
Aipij = (fii — tpj and Atp ik = ipi — ip k . The cumulant in terms of these relative angles 
can be obtained by summing Eq. 1 over all phase space. 

C3(Atpij, Aipn.) = J C3(ipi,(pj,tp k )5(Atfij — ipi + ipj)5(Aip ik - ipi + ipk)<bpidtpjd<pk 

The cumulants must be corrected for finite detection efficiencies, and averaged over 
collision centralities in heavy ion collisions. While this can be accomplished in a variety 
of ways, we advocate a simple technique based on the hypothesis that the efficiencies 
for simultaneously detecting multiple particles can be factorized as a product of the 
efficiencies for detecting each of the particles. Labeling the uncorrected cumulant as 
"raw", the efficiency corrected cumulant may be obtained as follows: 



C(A V A V ) - p i a _l {i) p ^ (j) P '"!l (fc) CI™' ' Ay ' fc) 
3 ^ e(j) e(k) pipip ^ ( Alpij , A Vik ) 

where the overbar indicate averaging over all azimuths, and Pl p 1 p[ J aw ^ (A<pij, A<pik) 
is calculated as follows: 

(raw) , A a \ r (raw) , \ (raw), (raw), * 

P1P1P1 {Apij,A(p ik ) = J p\ (VilPi Wj)P\ (f>k) 
x6(Aipij ~<fi + <pj)5(A<p ik -<Pi + <p k )d<fidipjd(p k 

Fig. f shows examples of three-particle density, combinatorial terms and 3-cumulants 
for 10-30% Au + Au collisions. 



3. Flow, Jets, and Conical Emission Discrimination 

Processes such as jet production, deflected jet production, and conical emission lead 
to distinct kinematical signature in C${Aipij , Aip ik ). In E], we modeled jet production 
within a Gaussian approximation, and found the away-side jet 3-cumulant component is 
proportional to exp ^— (A<p\2 — n) 2 /2a 2 2 — (Aipvi — it) 2 /2(t^ 3 — (Aip23) 2 /^^s) with 
an amplitude determined by the number of jets per event, and the number of associated 
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Fig. 1. Examples of (a) 3-particle density, (b,c) combinatorial terms, and (d) 3-cumulant obtained 
for Au + Au collisions at RHIC. Angles expressed in degrees. See text for definitions. 



particles per jet. The widths 012 and (T13, somewhat larger than 023 for in vacuum 
jets, are expected to significantly increase for in medium deflected jets while CT23 which 
measures the "intrinsic" width of the jet changes only modestly. Mach cone emission 
is, in principle, easily distinguishable from jet production as it leads to four jacobian 
peaks in C3 at A(pij,Aipn, = tt±8m,'k±8m where 8m is determined by the sound 
velocity in the produced medium . The shape and width of the peaks shall depend on 
the details of the production process and may be influenced by radial flow effects. Fig. 2 
illustrates the sensitivity of the technique based on a simulation involving near side jets 
with a width of 0.2 radian and conical emission at 8m with 0.4 radian width. It shows 
the 3-cumulants and projections along the main and alternate diagonals in slices of 40 
degrees. Plots were generated with a sample of 8 million events. Particle production was 
carried out with Poisson generators such that the event multiplicity ranges from 300 to 
600, with an average of 0.23 jets per event, and averaged associated multiplicities of one 
for the jet tag particle "1" and two for associates "2" and "3", in the near side jet and 
cone. 

The observation of conical emission and study of jet structure is complicated by 
the presence of irreducible flow components in the 3-cumulant. Assuming "background" 
particle emission is correlated to the reaction plane orientation and can be described 
with a Fourier decomposition e.g. rx 1 + 2 v n c °s (n (ipi — tp)) where v n and tp are 

n 

respectively flow coefficients and reaction plane angle, one finds the 3-cumulant shall 
contain contributions involving "off-diagonal" terms of the form: 



^ v p {i)v m (j)v n (k) X 



.„ cos (ptpi - mipj - ntpk) 
tP+n cos (— pifii + rmpj - n< 



+5 

_ +$n,rn + k COS (~p<Pi - mifj 



n<fk) 
ntpk) 



(5) 



Based on existing measurements of flow coefficients (see for instance HH), we expect such 
contributions to be dominated by terms of order V2V2V4 as indeed found in Fig. 1 (See 
El for details). 

The extraction of conical emission signal may be further complicated by correlated 
jet emission with the reaction plane. Indeed, one expects jet quenching or attenuation 
shall depend on parton path length in the medium. Given the finite spatial eccentricity of 
the medium produced in mid-central collisions, this implies jet emission can be correlated 
with the reaction plane orientation. Again using a Gaussian formulation for the jet 
profile, one finds that the 3-cumulant should contain jet-jet and jet-background terms 
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of the form: 

V2(jet)v2(bckg)exp(-Aipj j /2a? j ) cos (2(a4ipi + ajtpj - <p h )) (6) 

where the coefficients a; and a,j are determined by the relative widths of the jets. Fig. 3 
illustrates the shape of 3-cumulant resulting from such terms where we used = aj = 
0.5. Clearly, the reaction plane correlation produces an away-side narrow modulation, 
which should not be confused with three-particle jet contributions given its characteristic 
cosine dependence. 




Fig. 2. Top: Illustration of conical signal expected 
based on monte carlos simulation described in the 
text. Middle/Bottom: 40 degree slice projections 
along the main and alternate diagonals. 



Fig. 3. Top: Jet-Flow 3-cumulant component ob- 
tained with Eq. 6 for = aj = 0.5 plotted in ar- 
bitrary units. Middle/Bottom: Projections along 
the main and alternate diagonals. 



4. Summary 

We presented a new technique based on three-particle cumulants for the study of jet 
structure and searches for conical emission in heavy ion collisions at RHIC and LHC. 
While the cumulant technique was introduced for three-particle azimuthal correlations, 
it can be trivially extended to include correlations in rapidity space as well. We showed, 
with simple models of jet, flow, and conical emission, that it should be possible to 
distinguish these processes using three-particle correlations. Additional effects due to 
momentum conservation discussed inEl s hould also have distinct features in three-particle 
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cumulants. We thus believe that searches for conical emission with this type of observable 
should yield a definite and relatively unambiguous answer after proper treatment of flow, 
jet-flow background correlations, and momentum conservation effects. 
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